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I. INTRODUCTION 



It is widely accepted that low energy gluons are re- 
sponsible for the distinctive features of QCD such as con- 
finement and chiral symmetry breaking. Due to strong 
interactions, precision calculations of gluonic excitations 
are sparse and experimental evidence for gluonic excita- 
tions is still murky. This is because dynamical quarks 
and open channels need to be considered before compar- 
ing with experimental data. For the purpose of getting 
insights into the non-perturbative features of QCD it is 
however possible to design idealized systems that are sen- 
sitive to the pure Yang-Mills sector and to bypass many 
of the complications of the full QCD. A gluelump is an ex- 
ample of such an idealized state. It is a state of the gluon 
field bound to a static, localized octet source which can 
be constructed, for example as the quark-antiquark pair 
placed at an origin. Since the quark and the antiquark 
are at zero relative separation the system is rotationally 
symmetric, invariant under parity and charge conjuga- 
tion. Thus gluelump states can be classified by the same 
J PC quantum numbers as ordinary mesons. The spec- 
trum of gluelumps have been obtained from lattice gauge 
simulations and shows an unusual ordering of the various 
J PC levels [1, 2]. The ground state has J PC = l+~ and 
the first excited state has J PC = 1~~ . If the gluelump is 
to be interpreted as a state of a single (effective, quasi) 
gluon bound to a static source, then one would expect 
the 1 state to be below 1^ , since the former corre- 
sponds to a single gluon in an s-wave orbital and the 
latter to a gluon in a p-wave orbital. The following two 
levels have J PC = 2 and 2 H respectively and simi- 
larly correspond to a reverse ordering of the p-wave or- 
bital (2 H ) and the d-wave orbital (2~~). It is clear that 
a simple constituent model of the gluon bound to the 
QQ with a central potential [3] will not reproduce the 
lattice spectrum, since the centrifugal barrier orders the 
spectrum according to the orbital angular momentum. 
Another popular model for low lying gluonic excitations 
is based on the strong-coupling limit and considers glu- 



ons as excitations of the chromoelectric flux tube [4] . To 
the best of our knowledge the model has not been used 
to compute the gluelump spectrum but it does repro- 
duce the ordering of the few lowest levels of the excited 
adiabatic potentials between the QQ sources [5-8]. The 
adiabatic potentials give gluon energies for finite sepa- 
ration between the quark and the antiquark and reduce 
(modulo a zero point energy) to the gluelump spectrum 
in the limit of vanishing QQ separation. The correct or- 
dering of the adiabatic potential in the flux tube model 
is however somewhat artificial as it relies on a postu- 
lated intrinsic charge conjugation to be assigned to the 
flux tube. Furthermore, while there is a clear evidence 
that at large separation between the QQ sources the flux 
tube develops [9] and gluonic excitations are effectively 
those of the Nambu-Goto string [6], on the contrary it 
is not clear if the flux tube picture is correct when the 
QQ separation is taken to zero. In fact, the splitting 
of adiabatic potentials from lattice simulations for small 
QQ separations may be more constituent-like rather than 
string-like [10]. The gluelump spectrum may therefore be 
used to discriminate between alternative pictures of the 
gluonic excitations far from the adiabatic limit. Another 
reason for why it is important to pin down the charac- 
teristics of gluonic excitations for QQ separation smaller 
than ~ 1 fm is because of their relevance to the dynamics 
of heavy hybrid quarkonia. These are physical states and 
there already might be experimental evidence for such 
states [11]. Finally mixing between the QQ and the hy- 
brid component of a heavy quarkonium is responsible for 
relativistic corrections to the static QQ potential which 
is relevant to normal quarkonium spectroscopy [12-15]. 

In this paper we examine the gluelump spectrum in the 
variational approach based on the Coulomb gauge QCD. 
A variational vacuum wave functional constrains single 
gluon properties and in particular leads to an effective 
gluon mass. This enables truncation of the gluon Fock 
space and makes it possible to identify low lying gluelump 
states with the bound states of a single quasi-gluon. In 
Section II we present the details of the Coulomb gauge 
description of the gluelump spectrum and give the nu- 
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merical results. Summary and outlook are summarized 
in Section III. 



II. 



THE GLUELUMP SPECTRUM 



The variational approach to Coulomb gauge QCD [17] 
has recently been studied in [18-25]. Here we will 
only briefly summarize the method and give the rele- 
vant equations for the study of gluclump spectrum. The 
Coulomb gauge eliminates unphysical degrees of freedom, 
provided the domain of transverse field variables, A a (x), 
(V • A a (x) = 0) is restricted to the fundamental modular 
region (FMR) . The fundamental modular region is a sub- 
set of the Gribov region [19, 26] defined as the set of fields 
for which the determinant of the Faddeev-Poppov (FP) 
operator is positive and the FMR is a set of unique field 
representatives lying on a gauge orbit and satisfying the 
coulomb gauge condition. Restriction to the fundamen- 
tal modular region is hard to implement, however it was 
argued in [27] that the bulk contribution to the func- 
tional integrals over the field variables comes from the 
common boundary of the FMR and the Gribov region, 
and restriction to the Gribov region can be implemented 
within the variational approach. Within this approach 
the vacuum wave functional is chosen to be of the form 



*[A] = Ncxp (-± J dxA a (xMx-y)A a (y)) , (1 



with the gap function to determined by minimizing the 
expectation value of the Coulomb gauge Hamiltonian 
H = H[A,-iS/8A] 
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Here J = Dot [— V • D[A]] is the determinant of the 
Faddeev-Popov operator; D = D Q f, = <5 a bV + gf ac bA c 
is the covariant derivative in the adjoint representa- 
tion. The Coulomb gauge Hamiltonian H contains di- 
rect interaction between color charge densities, p a (x) = 
^(x)T a V(x) + /a hc n h (x) • A c (x), 



He 



J dxdyp a (x)K(x,a;y,b)p b (y) 



(3) 



where 



K(x, a;y, b) = {[V • Df 1 (-V 2 ) [V • Df 1 } (4) 

The inverse of the two covariant derivatives, when ex- 
panded in power series in the coupling constant, g, pro- 
duces an (infinite) series of interactions that couple an ar- 
bitrary number of transverse gluons to the quark and/or 
gluon sources via the bare Coulomb potential — 1/V 2 . 
When these gluons are integrated over the vacuum wave 
functional of Eq. (1), they dress the bare Coulomb po- 
tential and it results in an effective interaction which is 



approximately linear for large separations between the 
sources [21, 22, 24]. This effective Coulomb potential is 
defined as 



V C L(p)6ab = - / dxe ipx (K{x,a;0,b)) 



(5) 



with 

(K(x, a; 0, b)) = J L>AJV kx A"(x, a; 0, fe)|*[A]| 2 . (6) 

The subscript CL means that Vcl is expected to have 
a short distance, Coulomb-like component and a long- 
distance, linear part. The Dyson series for the effective 
interaction, Vcl is divergent in the ultraviolet and the 
bare coupling is rcnormalizcd by fixing the slope of the 
potential to agree with the lattice data for the Wilson 
loop. There is, however, a difference between the renor- 
malization group equations for the coupling obtained in 
the variational approach and that of the QCD beta- 
function. This is to be expected since the variational 
approximation includes only two-body correlations ( just 
like BCS) and, for example does not take into account 
propagation of transverse gluons [18, 22]. Furthermore, 
even for exact vacuum the expectation value of He is 
not the same at the energy of the QQ state [28]. This is 
because the state of a QQ pair added to the vacuum is 
not the same as the true eigenstate of the QCD Hamil- 
tonian in presence of the QQ pair. All this lead to an 
ambiguity in fixing the long range part of the QQ po- 
tential. Fortunately we expect that this will not play an 
important role for small, and in the case of gluelumps, 
vanishing QQ separations, and in general for low energy 
gluonic excitations that have wave functions confined to 
a region of ~ 1 — 2 fm that are not too sensitive to the 
asymptotic, long range behavior of the confining poten- 
tial. Finally, ambiguities in the renormalization group 
equations of the variational methods, have so far pre- 
vented from obtaining a unique solution for uj by using 
the gap equation [23, 29, 30]. In particular it is expected 
that at low momentum, oj(p), (the Fourier transform of 
lu(x — y)) has momentum dependence identical to that 
of the expectation value of the curvature, xip) defined 
through the determinant of the Faddeev-Popov operator, 
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v <SA a (x)A h (0)' 



(7) 



while at large momenta, oj(p) ~ p. There are indications 
from lattice computations that to(p) is different from zero 
as p — > but the lattice results cannot resolve between 
Lu(p — > 0) finite and u)(p — > 0) — > oo [31, 32]. In our 
analysis we will thus study how the gluelump spectrum 
depends on the gluon gap function uj. 

The gap function to is not the same as the single quasi- 
gluon energy, E(p), although they are closely related. 
For uj(p) satisfying the gap equation, Eq. (2), the single 
particle energy is given by 



E{p) = u(p) + £(p) 



(8) 



3 



with the self-energy given by 



Vb£(|p-k|) l + (p-k 



>(p) 



ui(k)' 
(9) 

The long range nature of the Coulomb kernel K, Eq. (4), 
leads to strong IR enhancement in Vcl(p), Vcl{p — > 
0) — > oo, which in turn makes the integral in Eq. (9) 
divergent for p ~ k. This is a signature of color confine- 
ment in the mean field which leads to infinite energies for 
colored states. The total energy of a color singlet state, 
in particular gluclumps, is however finite, as will be seen 
below, with all theIR divergences canceling among vari- 
ous contributions. 



A. Fitting the potential 

The solution of the Dyson equation for the effective 
interaction Eq. (5) can be written in terms of the square 
of the expectation value of the inverse Faddeev-Popov 
operator, d(p), 



d(p) 

p2 



Sab 



dxe ipx ( 



V D 



)(x,a;0,6) 



(10) 



and the Coulomb form factor, f(p) which measures the 
difference between the square of the expectation value of 
the FP operator and the expectation value of the square 
of the FP operator which defines the Coulomb kernel 
Eq. (4), 



Vcl(p) 



d 2 (p)f(p) 



(11) 



As discussed above, the mean field approximation, by 
neglecting transverse gluon loops, leads to a high mo- 
mentum behavior of Vcl{p), which is different from that 
of QCD. For p -> oo V CL (p) oc l/[p 2 log 3/2 (p)]. In the 
low momentum regime, to reproduce the linear confine- 
ment, Vchip) should grow as 1/p when p — > 0. Unfor- 
tunately, in the mean field approximation common so- 
lutions of the Dyson equations for d, f and to do not 
quite reproduce such a strong IR enhancement of Vol- 
In the so called angular approximation [22] it is possible 
to obtain a solution with both d(p) and f(p) divergent 
in the limit p — > resulting in Vcl(p) oc 1/p 15 / 4 . It was 
found in [25, 29], however, that without the angular ap- 
proximation the three Dyson equations for d, f and u> 
do not have solutions such that both d(p) and f(p) are 
divergent in the IR (as p — > 0). The IR behavior of the 
solutions is controlled by the value of the renormalized 
coupling g — > g{p). In the mean filed approximation it 
is found that there is a critical coupling g c {p) for which 
d(p) is IR divergent, while for all g(fi) < gdn) both d(p) 
and f(p) saturate in the IR. These sub-critical solutions 
do not, in a strict sense lead to confinement, however, 
for all practical application Vcl{t) (defined as a three- 
dimension Fourier transform of Vcl(p)) computed from 



these sub-critical d and / is very close to a linear function 
of r for a r < 10 fm. We will thus approximate Vcl{t) by 
a linear potential for r in this range plus a short range, 
Coulomb piece modified by the loop corrections as dis- 
cussed above. In momentum space, 



V CL (p) = V c (p) + V L (p) 



(12) 



with the Coulomb part, Vc, and the linear part, Vl, given 
by 



V c (p) = 
Vl(p) = 



4ira(p) 

87T& 



*(P) 



4ttZ 



\ 1V QCD 



(13) 



The Fourier transform of Vl which defines the coordi- 
nate space potential is ill-defined due to the strong IR 
singularity in Vl , and leads to an undetermined constant 
in Vcl{ t ) = Vo(r) + Vt(r). In computations of Vcl{t), 
including lattice, this constant is removed by computing 
a difference Vcl{t) — Vcl{tq) at some fixed tq: 



V L (r)-V L {r ) 



dp 

(2^js 



8tt6 



(14) 



If we were to remove the logarithm in the Coulomb part 
a(p) — * a = const, then also the Coulomb part would 
have a simple position space representation 



V c (r) 



dp 



Aira 

P 



a 



(15) 



In the following, however, we will keep the full expression 
for a(p) and compute the position space potential Vc(r) 
numerically 



V c (r) = -(4tt) 



p 2 dp Aitaip) 
—^Jo(rp) — 



(27T) 



P- 



(16) 



We fix the four parameters, Z, c, Kqcd, b, by fitting the 
position space potential Vcl(t) — Vcl{to) to the lattice 
data [33] with r - 1/450 MeV. The result of the fit 
is shown in Fig. 1. Unfortunately the low-r lattice data 
does not go to small enough distances to fix the Coulomb 
potential uniquely. Thus in the fit we fixed Aqcd = 
250 MeV, and then obtain 



b = 0.204 GeV 2 , Z = 5.94, c = 40.68. 



(17) 



B. Basis of Gluelump states 

The mean field vacuum leads to a basis of quasi- 
particle states which can be used to express any state 
of the Yang-Mills field. The quasi-particle creation and 
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r = R/r„ 



FIG. 1: Lattice data and fit results. The errors on the lattice 
data are taken to be Se = 0.05 and the fit gives x /P-d-f = 
1.85 

annihilation operators are defined in terms of the field 
variables 

a-w = £/^7^<>a,«)e(M) 

+ a t (-k J A,a)e*(-k 1 A)]e aMC 
n«(x) = -i£ / ^ ?V /^[a(k l A,«) 6 (k 1 A) 
- a t (-k,A,a)e*(-k,A)]e ik - x , 

(18) 

where II is the variable canonically conjugated to A (in 
the Schrodinger representation used in Sec.ll II a (x) = 
— iS/SA a (x)) and e(k, A) is the wave function of a state 
with hclicity-A, A = ±1. With w satisfying the gap equa- 
tion, the Hamiltonian rewritten in terms of the particle 
operators a, has a diagonal one-body part, i.e it con- 
tains only the ladder operator 

K = E / (^3 E (*) at ( k - A ' «)«( k ' A ' °). ( 19 ) 

while operators of the form aa and a^aJ are removed by 
the gap equation. The single particle energy E{k) is given 
by Eq. (8). Since (after IR regularization) the quasi- 
particles are massive, it is reasonable to expect that light 
gluclump states will contain a minimal number, which is 
one, of quasi-particles constrained by various symmetries. 
The state of a single-quasigluon gluclump is given by 

\ jPM ) = E/ ^aQl^3«+(k,A,a)|0) 

x -Jpl—^kA). (20) 

The SU(Nc) generators, T a couples the iVc-quasi-gluons 
to the quark-antiquark state in the adjoint representation 



resulting in a color-singlct gluelump state (Cf = (Nq — 
\)/2Nc)- Quark spin and position quantum numbers are 

implicit in the quark (Qt) and antiquark (Q^) creation 
operators. The orbital wave function describes spin-1 
quasi-gluon with helicity A = ±1 projected onto a state 
with total spin- J, projection-M and parity-P and it is 
given by 

< ( kA) = J^±± J2 X&DtiMV P (k). (21) 

A'=±l 

Here xw distinguishes between the TM and TE modes, 
XXX' = Sxx'/VZ for TM (natural parity) and xaa' = 
\6\\>/V2 for TE (unnatural parity) modes, respectively, 
and ip J (k) is the orbital wave function which has to be 
determined by diagonalizing the QCD Coulomb gauge 
Hamiltonian in this single quasi-gluon subspace. It fol- 
lows from the properties of the Z?-functions that the ra- 
dial wave functions are normalized by 

/ ^3<P jP (W J ' P 'W= S JJ' S r?'- ( 22 ) 

Finally we note that all gluelump states built from single 
quasi-gluon have negative charge conjugation. 

C. Hamiltonian matrix elements 

Having constructed the basis which complies with the 
global symmetries of the Hamiltonian, we can compute 
the matrix elements in the single quasi-gluon subspace. 
Digitalization of the Hamiltonian matrix leads to a one- 
dimensional Schrodinger equation for the radial wave 
functions, ip J (k). 

E JP i> jP (p) = [Lu(p) + X(p) + 2™ Q + X QQ + V s ]t JP (p) 
+ J^V jP ( P ,k)t JP (k) (23) 

Here ui + E is the quasi-gluon kinetic energy given in 
Eq. (9), Eqq is the sum of source (quark an antiquark) 
self energies, 

*qQ = -Cp J J^sVcUD, (24) 

and V s is the potential energy of the octet QQ pair at 
zero separation 

v °'-iklw? v ° L(l) - (25) 

The QQ self energy and the V s are infinite. The in- 
finities come both from the short distance part of the 
Vc potential and from the long range part of the linear 
potential Vl- Part of the UV divergence should be re- 
moved by the quark masses (toq) and part needs to be 
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subtracted because the QQ pair at zero separation has 
infinite chromo-electrostatic energy. The IR divergence 
is a manifestation of the long-range nature of the inter- 
action and as long as the system is color-neutral, which 
is the case here, should cancel with similar contributions 
from gluon-sclf energies and residual interactions. In or- 
der to remove the infinite chromo-electrostatic energy the 
gluclump energy is defined as [2] 



Eq = E Jp - E Q q -Vc + Vc, 



(26) 



where Eqq is the total energy of the color-singlet, static 
QQ source 



with 



E QQ = 2m c 



V" = C F 



Sqq + V°, 



dl 



V CL (l) 



(27) 



(28) 



being the potential energy of the singlet QQ pair; Vq and 
Vq are the short distance parts of the octet and singlet 
QQ potential energy, 



1 



c 



d\ 



2N C J (2tt) 3 



Veil) 



(2tt) 3 



(29) 



respectively. Combining these in the Schodingcr equation 
for the gluelump energy we obtain 



The modified gluon self energy £' is now finite in both 
IR and UV and, as will be seen below, AV removes the 
IR divergence from the matrix elements of the potential 
between the quasi-gluon and the quark or the antiquark. 
The potential V J is a sum of two- and three-body in- 
teractions. 



V 



(33) 



The two body interaction acts between the quasi-gluon 
and the quark or the antiquark while the three-body 
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FIG. 2: Two -(left) and three- (right) body potential between 
the quasi-gluon and the QQ pair. The static quark and the an- 
tiquark are represented by blobs. The dashed line represents 
the effective Vcl potential (left) or the effective three-body 
interaction from Eq. (37) (right). 



E& t' (p) = Kp) + £»]V/ (p) 
dkk 2 , rJ t 



where 



S'(p) 



and 



Nc 
2 



(2tt 



dk 



|3 



V J ' (p, k) + AV J (p, k) V' 7 (k) 

(30) 



V CL (\p-k\) 



1 + (p ■ k) 2 to(p) 
2 Lu(k) 



- 1 



AV J = - 



Nc { ^S(k-p) 



dl 



(31) 



■V(l). (32) 



terms involve all three particles. Both terms originate 
from the Coulomb kernel of the QCD Hamiltonian and 
are depicted in Fig. 2. In the two body interaction all 
gluons emerging from the Coulomb kernel, Eq. (4), are 
integrated over the vacuum wave functional, producing 
the effective interaction Vcl which acts as a potential 
between the quasi gluon and the quark or the antiquark. 
In the three-body term V^ , the gluon in the gluelump 
is contracted with one of the gluons in the Coulomb ker- 
nel, producing an effective interaction which involves all 
the three constituents, the gluon, the quark and the anti- 
quark. The two-body interaction, combined with AV J , 
is given by 



Vf{p,k) = + 



N c 



M*0 , MP) 



2/4-1 r 
4tt£ [1 - P(-l)<] gj^H, /0| Jl) 2 J d{p ■ k)V CL (k p)P l (p ■ k) (34) 



r 



where 

V^ L (|k-p|) = V Ci (|k-p|)-(2^) 3 ( S(k-p) 



is IR finite. For natural parity states, P = {— l)' 7 , the 
ji two-body potential reduces to 

Vl(1) 



O) 



(35) 
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V 2 J (p,k) = + 



M*0 Mp) 
u(p) + V 



4tt 



2J+1 



d(p ■ k)V£ L (k - P )[(J + • k) + JP J+1 (p • k)], 



and for unnatural parity states, P = ( — 1) J+1 , 



V? 



ip,k) = +— 




The three-body interaction is by itself IR finite and its 
matrix elements are given by [34] 



J 



An / d(£.k)V& £ (k-p)Pj(p.k). 



V(p,q) 



N% (47T) 2 
4 ^/MpjMqj 



X) [i--^- 1 ) 1 *] (i,o > ii|Ji)(jpO,ii|Ji)(z a o,io|jo>(i p o,io|jo) 



(2Z g + 1)(2Z P + 1) / dfcfc 



(2 J +1) 2 



(2tt) 3 



d(q • k)d(p • k)fc 2 X(|q + k|, |k|, |k + p\)P lq (q • k)P ip (p • k) 



r 



(36) 



where the three-body kernel is 

x [/(p)d(p) + /(fc)d(fc) + /(?)<*(?)] . (37) 

Each of the three terms in square bracket is given by a 
product of a Vcl potential, Eq. (11), in either p, or k 
or q multiplied by bare Coulomb potentials in the other 
two momenta corrected by the ghost form factors. In 
numerical computations we approximate these two mod- 
ified coulomb potentials by Vc times a scale factor, n, 
which is fitted to the data. Thus, finally 



K(p,q, k) 



K{V C L(p)Vc(q)Vc(k) 
V C (p)V C (q)V CL (k)}. 



Vc(p)V C L(q)Vc(k) 



D. Numerical results 



Since the quasi-gluons have hclicity ±1 the Wigner-D 
functions in Eq. (20) restrict the spin of the low lying 
gluelumps, saturated by a single quasi-gluon, to J > 1. 
This is in agreement with lattice computations [1, 2]. 
The two lowest J-states, arc J p = 1~ and J p = 1 + . 
For the natural parity the gluclump state is a mixture of 



the I = J — 1 = 0, s-wave, and I = J + 1 = 2, d-wave, 
quasi-gluon state, which is explicit in the form of the two- 
body interaction in Eq. (36). The unnatural parity state 
contains the quasi-gluon in the I = J = 1, p-wave state, 
cf. Eq. (36). It is therefore expected that the two-body 
interaction alone would lead to the energy of the natural 
parity, J p — 1~, state to be lower than that of the unnat- 
ural parity, J p = 1 + , state. The lattice computations, 
however find the J p = 1 + state to have lower energy 
than the J p = 1~ state. This is due to the three-body 
potential . Indeed, for J = 1 and positive parity, in 
Eq. (36) the [1 — P(— l) lq ] term implies l q = odd and the 
Clebsch-Gordon (CG) coefficient (Z g 0,10|J = 1,0) van- 
ishes. For J = 1 and negative parity, the product of all 

the four CG coefficients is positive. Thus the three-body 
(38) 

interaction does not contribute to the J p = 1 + state, 
while it increases the energy of the J p — 1~ state. The 
three-body interaction can therefore change the ordering 
of the 1~ and 1 + state and push the 1" state above the 
1 + state. Similar inversion can take place for higher J 
states. In general V% vanishes for all unnatural parity 
states. Thus, while for given J the two-body potential is 
expected to produce the natural parity state with lower 
energy than the unnatural parity partner, the three-body 
interaction will add to the energy of the natural parity 
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FIG. 3: Three models for the gap function w(p). 



state, possibly pushing it above the unnatural parity one. 
This happens in the lattice data at least for the first five 
states. In fact 1 + is below 1~. The next state is the 2~ 
and is below the 2 + state, which is followed by the 3 + . 
The sixth state is the ++ which, having both positive 
charge conjugation and J = 0, in the quasi-gluon picture 
would require two gluon constituents. 

For each J p gluclump state the Schrodinger equation, 
Eq. (30), is solved by discretizing the momentum coordi- 
nate through a Gauss-Legendre grid. We examine sensi- 
tivity of the computed spectrum to the form of the gap 
function lu. In particular we first try two cases. One with 
Lu(k — ► 0) = finite and the other with u>(k — ► 0) — * oo. 
The former has been found as an approximate solution to 
the Dyson equations in Ref. [22, 23] and is referred to as 
modcl-1. The latter has been advocated in Ref. [29, 30] 
as the solution of the IR limit of the Dyson equations 
which include the curvature but set the Coulomb form 
factor f(p) = 1 and is referred to as model-2. These gap 
functions are shown in Fig. 3. The gluelump spectrum 
for these two gap functions is displayed in Fig. 4 together 
with the lattice results from Ref. [2]. The vacuum wave 
functional built from these gap functions systematically 
overestimates the masses of the gluclumps by approxi- 
mately 0.7GeV. The inversion of the levels, due to the 
three-body interaction, is clearly visible with 1 _ above 
1 + and 2~ near but still above 2 + . In the next step we 
varied u in the attempt to bring the 1 + state to be in 
agreement with lattice mass. We have found that the over 
mass scale is, as expected, mostly sensitive to the value 
of uj{p) in the region 0.5 GcV < p < 1.0 GcV. We thus 
chose a new functional form for the gap function which 
reduced uj in this momentum region. This is referred to 
model-3 in Fig. 3. The corresponding gluelump spectrum 
is shown in Fig. 5 where we show the result of the calcu- 
lation without the three-body interaction (dashed lines) 
and the full spectrum (solid lines). As discussed above 
the two-body potential leads to the inverted spectrum 
with 1~ below 1 + and 2 + below 2 _ while the three-body 



FIG. 4: Gluelump spectrum from Eq. (30) with uj from models 
1 and 2 shown in Fig. 3 and n = 0.6 (Eq. (38)) compared with 
results of lattice computations [2] . 



E a [GOV] 



lattice 

■ ■ model 3 two-body interaction only 
— model 3 two and three-body interactions 



FIG. 5: Gluelump spectrum with u>(p) from model-3, shown 
in Fig. 3, and k — 0.4 (Eq. (38)). The dashed lines correspond 
to the spectrum obtained with two-body interaction only and 
the solid lines correspond to the full calculation. 



interaction is operative for the natural parity states only, 
bringing the Coulomb gauge spectrum in agreement with 
the lattice spectrum. 



III. SUMMARY AND OUTLOOK 

The quenched gluclump state is the simplest system 
that can be used to define a physical state of a gluon. 
Within the framework of Coulomb-gauge QCD with a 
non-perturbative ansatz for the vacuum wave functional 
gluons correspond to quasi-particle excitations and low 
lying gluclumps are states of a single quasi-gluon whose 
color is neutralized by an external static source. The 
possible quantum numbers of such states match with 
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those found in lattice computations. Furthermore, order- 
ing of these levels is reproduced, and it turns out to be 
non-trivially related to the non-abelian structure of the 
Coulomb interaction in the Coulomb-gauge which leads 
to new (three-body) interaction in the Schrodinger equa- 
tion for the gluelump energies. The shape of the gap 
function in the low (jp < 1.0 GeV) is primary relevant 
for setting the energy scale of the gluclumps. Because of 
the ambiguities in the rcnormalization program of vari- 
ational Coulomb gauge calculations originating from the 
truncation of the Fock space, these findings can be used 
to further constraint the Dyson equations of the mean 
field Coulomb gauge calculations. 
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